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Abstract: 

We provide a distributed coordinated approach to the stability analysis and control design of large- 
scale nonlinear dynamical systems by using a vector Lyapunov functions approach. In this formulation 
the large-scale system is decomposed into a network of interacting subsystems and the stability of the 
system is analyzed through a comparison system. However finding such comparison system is not trivial. 
In this work, we propose a sum-of-squares based completely decentralized approach for computing the 
comparison systems for networks of nonlinear systems. Moreover, based on the comparison systems, we 
introduce a distributed optimal control strategy in which the individual subsystems (agents) coordinate 
with their immediate neighbors to design local control policies that can exponentially stabilize the full 
system under initial disturbances. We illustrate the control algorithm on a network of interacting Van der 
Pol systems. 
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1. INTRODUCTION 

Distributed coordinated control has recently provided power¬ 
ful control solutions when the conventional centralized meth¬ 
ods fail due to inevitable communication constraints and lim¬ 
ited computational capabilities. Paradigmatic examples are pro¬ 
vided by cooperative and coordinated control for autonomous 
multi-agent systems (see Bullo et al. (2009)) or large scale inter¬ 
connected systems (see Zecevic and Siljak (2010)). Distributed 
coordinated control uses local communications between agents 
to achieve global objectives that reflect the desired behavior of 
the multi-agent system. Usually, a two-level hierarchical multi¬ 
agent system is employed, which consists of upper level agent 
for implementing coordinated control and lower level agents 
for implementing decentralized control. In this paper, we pro¬ 
pose to use this conceptual framework to design distributed 
coordinated control of large scale interconnected system us¬ 
ing vector Lyapunov functions (see Bellman (1962); Bailey 
(1966)) and comparison principles (see Brauer (1961); Beck- 
enbach and Bellman (1961)). The formulations using vector 
Lyapunov functions are computationally very attractive because 
of their parallel structure and scalability. However computing 
these comparison equations, for a given interconnected system, 
still remained a challenge. In this work we use sum-of-squares 
(SOS) methods to study the stability of an interconnected sys¬ 
tem by computing the vector Lyapunov functions as well as the 
comparison equations. While this approach is applicable to any 
generic dynamical system, we choose a randomly generated 
network of modified * 1 Van der Pol oscillators for illustration. 

* This work was supported by the U.S. Department of Energy through the 
LANL/LDRD Program. 

1 We choose the Van der Pol ‘oscillator’ parameters in such a way that these 

have a stable equilibrium at origin. 


This network is decomposed into many interacting subsystems 
and each subsystem parameters are chosen so that individually 
each subsystem is stable, when the disturbances from neighbors 
are zero. SOS based expanding interior algorithm (see Jarvis- 
Wloszek (2003); Anghel et al. (2013)) is used to obtain estimate 
of region of attraction as sub-level sets of polynomial Lyapunov 
functions for each such subsystem. Finally SOS optimization 
is used to compute the stabilizing control policies, based on 
linear comparison systems, such that the closed-loop network 
is exponentially stable under initial disturbances. 

Following some brief background in Section 2 we formulate 
the control design problem in Section 3. The sum-of-squares 
based distributed control algorithm is proposed in Section 4. In 
Section 5 we illustrate the control design on a network of Van 
der Pol systems, before concluding the article in Section 6. 

2. PRELIMINARIES 

2.1 Stability and Control of Nonlinear Systems 

Let us consider the dynamical systems of the form 

x{t) =f(x(t))+u t , t >0, /(0) = 0, (1) 

where x £ M" are the states, it, £ R" are the control input, 
/: R" —>• R" is locally Lipschitz and the origin is an equilibrium 
point 2 of the ‘free’ system, i.e. the system with no control 
(u, = 0). Let us first review the important concepts on stability 
of the equilibrium point of the ‘free’ system. 

Definition 1. The equilibrium point at the origin is called 
asymptotically stable in a domain fCl", 0gS>, if 
||x( 0)|| 2 £^ => lim ||x(f)|| 2 = 0, 

t—foo 

2 State variables can be shifted to move any equilibrium point to the origin. 






and it is exponentially stable if there exists b,c> 0 such that 

||x(0)|| 2 e^=>||x(0||2<ce“ 6, ||x(0)||2 Vf>0. 

Lyapunov’s first or direct method (see Lyapunov (1892); Slo- 
tine et al. (1991)) can give a sufficient condition of stability 
through the construction a certain positive definite function. 
Theorem 1. If there exists a domain ^£R", 0gS>, and a 
continuously differentiable positive definite function V:R" —>R, 
called the ‘Lyapunov function’ (LF), then the equilibrium point 
of the ‘free’ system at the origin is asymptotically stable if 
VV r /(x) is negative definite in 5L and is exponentially stable if 
VV f f(x) < —cV MxC3>, for some c > 0. 

When there exists such a V (x), the region of attraction (ROA) 
of the equilibrium point at the origin can be estimated as 

M:= {x€@\V(x)<l}, (2a) 

where, V (x) = V (x) /Y nax , and (2b) 

yno.x ar g max |x £ M" |y(x) < y} C (2c) 

For systems under some control action u t , the notion of ‘stabi- 
lizability’ becomes important. Specifically, we are interested in 
state-feedback control of the form u t = u, (x). 

Definition 2. The system (1) is called (exponentially) stabiliz- 
able if there exists a control policy u t = u t (x), t > 0, such that 
the origin of the closed-loop system is (exponentially) stable, in 
which case u t is called a (exponentially) stabilizing control. 

Courtesy to the works of Artstein (1983) and Sontag (1989), the 
concept of ‘control Lyapunov functions’ has been useful in the 
context of stabilizability. 

Definition 3. A continuously differentiable positive definite 
function V, : R" —> R is called a ‘control Lyapunov function’ 
(CLF) if for each x £ R”\{0}, there exists a control u t such that 

VV c 7 ’(/(x) + n,) <0. 

Similar definition holds for ‘exponentially stabilizing’ CLFs 
(see Ames et al. (2014); Zhang et al. (2009)). CLFs can easily 
accommodate ‘optimality’ in the control policies as well (see 
Freeman and Kokotovic (2008)). However, as with the LFs, it 
is often very difficult to find a CLF for a given system. 

2.2 Sum-of-Squares and Positivstellensatz Theorem 

In recent years, sum-of-squares (SOS) based optimization tech¬ 
niques have been successfully used in constructing LFs by 
restricting the search space to sum-of-squares polynomials (see 
Jarvis-Wloszek (2003); Parrilo (2000); Tan (2006); Anghel 
et al. (2013)). Let us denote by R [x] the ring of all polynomials 
in x £ R”. Then, 

Definition 4. A multivariate polynomial p £ R[x], x £ R", is 
called a sum-of-squares (SOS) if there exists £ R[x], i £ 
{1,... ,s}, for some finite s, such that p(x) = £T =1 hj(x). Fur¬ 
ther, the ring of all such SOS polynomials is denoted by E[x], 

Checking if p £ R[x] is an SOS is a semi-definite problem 
which can be solved with a MATLAB® toolbox SOSTOOLS 
(see Papachristodoulou et al. (2013); Papachristodoulou and 
Prajna (2005)) along with a semidefinite programming solver 
such as SeDuMi (see Sturm (1999)). SOS technique can be used 
to search for polynomial LFs, by translating the conditions in 
Theorem 1 to equivalent SOS conditions (see Jarvis-Wloszek 


(2003); Wloszek et al. (2005); Prajna et al. (2005)). An impor¬ 
tant result from algebraic geometry called Putinar’s Positivstel¬ 
lensatz theorem 3 (see Putinar (1993); Lasserre (2009)) helps in 
translating the SOS conditions into SOS feasibility problems. 
Theorem 2. Let JC= {x £ R" | k\ (x) > 0,..., k,„ (x) > 0} be a 
compact set, where kj £R[x], Vj £ {1 Suppose there 

exists a p £ joo + Y!j=iGjkj \ Go,Gj £ £[x] ,Vy j such that 
{x £ R"| /i(x) > 0} is compact. Then, if p(x) >0 Vx£JC, then 
P £ { Off +Lj Gjk j\ Go , Gj £ £ [x], V / } . 

In many cases, especially for the k,Vi used throughout this 
work, a fi satisfying the conditions in Theorem 2 is guaranteed 
to exist (see Lasserre (2009)), and need not be searched for. 

2.3 Linear Comparison Principle 

Before finishing this section, let us take a look at a nice result 
on the ordinary differential equations which helps form the 
framework of stability analysis of inter-connected systems via 
vector LFs. Noting that all the elements of the vector e At , t > 
0, where A = [aif £ R mxm , are non-negative if and only if 
ajj > 0. i Y j , the authors in Beckenbach and Bellman (1961); 
Bellman (1962) proposed the following result: 

Lemma 1. Let A = [a,y] £ R mxm have only non-negative off- 
diagonal elements, i.e. a l; >0, i / j. Then 

v(t) <Av(t ), t> 0, t' £ R", v(0) = vo, (3) 

implies v(t) < r(t), Vf > 0, where 

r(t) =Ar(t), t >0, r £ R", r(0) = v(0) = v 0 . (4) 

This result will henceforth be referred to as the ‘linear com¬ 
parison principle’ and the differential equation in (4) as the 
‘comparison equation’. 

3. PROBLEM DESCRIPTION 

The problem of interest for this work is to find state-feedback 
control u t = u, (x) that exponentially stabilizes a large nonlinear 
system (1). One approach could be to find a suitable CLF (Def¬ 
initions), using computational methods, e.g. SOS technique. 
However, as noted in Anderson and Papachristodoulou (2012), 
such an approach will quickly become intractable as the system 
size increases. Instead, we seek distributed stabilizing control 
policies by modeling the large dynamical system as an inter¬ 


connected network of m (> 2) interacting subsystems, 

Vi = 1,2,... ,m, 

: Xi = fi(xi) + u t j + gi(x), xi £R"', x£R" (5a) 
fi( 0)=0, (5b) 

gi(xi) = 0, Vx,- £ {x £ R"| xj=0,Vj^i} (5c) 

where, x = | { x f } , and n < ^"'_,n . (5d) 


We assume that the isolated ‘free’ subsystem dynamics /,- £ 
M[x,]"', and the neighbor interactions gj £ R[x]"' are vectors of 
polynomials. Further, u t j = u tJ (x,) is a time-dependent local 
state-feedback control policy, with each n,,- £ R[x,]"' Vf. It is 
assumed that the ‘free’ isolated subsystems as well as the 
‘free’ full system are (locally) stable. Note that, we allow over¬ 
lapping decomposition in which subsystems can have common 
state(s) Siljak (1978); Jocic and Siljak (1977). Let 

3 Refer to Lasserre (2009) for other versions of the Positivstellensatz theorem. 



JVi-.= {;'} U 


gi(x) Y 0 for some x 1 
with x k = 0 Vk^i,j J ’ 


(6a) 


and xr-=\J jeA; { x j} (6b) 

denote the set of indices of the subsystems in the neighborhood 
of (including the subsystem itself) and the states that belong 
to this neighborhood, respectively. 


The goal is to compute the distributed control m,,(x,) Vi so that 
the full interconnected system (5) is exponentially stabilizable. 


3.1 Comparison Equations and Exponential Stabiltiy 


Let us first review the stability of the ‘free’ interconnected 
system, i.e. when u t j = 0Vi. Stability of each of the ‘free’ 
isolated (i.e. zero neighbor interaction) subsystems 

Vie {1,2,...,m}, Xi = fi(xi),XiGM. n ‘. (7) 

can be characterized by computing a polynomial LF V)(x,-)Vi, 
and the corresponding estimate of the ROA as in (2). An 
SOS based expanding interior algorithm , (see Jarvis-Wloszek 
(2003); Anghel et al. (2013)), is used to iteratively enlarge the 
estimate of the ROA by finding a ‘better’ LF at each step of 
the algorithm. At the completion of this iterative algorithm, the 
stability of each ‘free’ isolated subsystem (7) is quantified by 
its LF Vj(xj), with a corresponding estimate of the ROA as 

:= {x t G R"' |Vj(xj) < 1}, Vf = 1,2,... ,m. (8) 

Let us further define the domain 

:= {x G R" | x, G Vi = 1,2,.. .,m } . (9) 

The equilibrium of the ‘free’ network at the origin corresponds 
to the zero level-sets, V,(0) = 0 Vi, and any initial condition 
away from this equilibrium would result in positive level-sets 
V/(x/(0)) = yf G (0,1] for some or all of the subsystems. 


An attractive and scalable approach for (exponential) stability 
analysis of the ‘free’ network uses a vector LF (see Bellman 
(1962); Bailey (1966)) 


V(x) := [Vi(xi) V 2 (x 2 ) ... V m (x m )] T (10) 

to construct a linear comparison equation (Lemma 1) whose 
states are the subsystem LFs (see Siljak (1972); Weissenberger 
(1973); Araki (1978)). The aim is to seek an A = [a tJ \ G R mxm 


and a domain C & 0 , such that 

V (*) | u t i= 0 Mi < AV (x), Vx G 9 C (11a) 

where, a\j > 0 V; Y j . (lib) 

A = [atj\ is Hurwitz, and (11c) 

is invariant under the dynamics (1), (lid) 


Wf {fi(xi)+gi{x)) 


and V ( x K^0Vi = 


(lie) 


L ( fm{Xm ) + gm(x)) J 

If there exist a ‘comparison matrix’ A = [a i; ] and 0 C A? 0 
satisfying (11), then any x(0) G would guarantee exponential 
convergence of V(x(t )) to the origin thereby implying exponen¬ 
tial convergence of the states themselves (see Siljak (1972)). 


3.2 Exponentially Stabilizing Control 


The comparison principle can be used to design distributed 
controllers m^, (x,)V/ that exponentially stabilize the nonlinear 
network (5). In Section 4, we propose an SOS based algorithmic 
approach in which each of the subsystems coordinates only 


with its immediate neighbors , j G ,yf/\{;}, to compute a 
local and ‘optimal’ stabilizing control «, ,■. 

We propose that the LFs for each ‘free’ (no control) and 
isolated (no interaction) subsystem (7) be pre-computed and 
communicated to the neighbors. Given any initial condition 


x(0) G^° we define the domain 

& := {x G | Vi( Xi ) < Vj(xi{ 0)) = Vi} . (12) 

Then any distributed control m,,(x ; ) Vi satisfying 

V(x) < AV (x), Vx G 9 C ^°, (13a) 

s.t., conditions (lib), (11c) and (lid), (13b) 

v (/1 (xi) + u ti i (xi) + gi (x)) 

where V(x)= : . (13c) 

VL m ( fm(Xm ) T U,. m (x m ) A gin (x)) 


is an exponentially stabilizing control policy. In addition to sat¬ 
isfying (13), the ‘optimality’ of the control could be ascertained 
by minimizing the applied control efforts. 

Remark Note that we do not explicitly compute a CLF (Def¬ 
inition 3), because of the computational burden in large-scale 
networks. Instead, we propose an algorithm to design stabiliz¬ 
ing control using the pre-computed subsystem LFs. 

4. DISTRIBUTED CONTROL ALGORITHM 

In designing the stabilizing control policies m / ( -V; in (13) the 
conditions (11c) and (lid) have to be satisfied, which essen¬ 
tially demands availability of network-level information. How¬ 
ever, the following two key observation can be useful in gener¬ 
ating equivalent subsystem-level conditions. 

Proposition 1. A matrix A = [a,y] G R mxm is Hurwitz if, for 
each * G { 1 , 2 ,.. . ,m}, au + 'Lj^i\ a ij\ < 0. 4 

Proof From the Gershgorin’s Circle theorem (see in Bell 
(1965); Gershgorin (1931)), for every eigenvalue A G C of the 
matrix A = [a i; ], 

3k G {1,2,... ,m} such that, |A — \ < \ ak i\ • 

Using Y.jjtk | a kj | < ~ a kki it follows that Re{A} < 0 . □ 

Additionally, we also note that (see Weissenberger (1973)), 
Proposition 2. The domain @ in (12) is invariant if Y!J= 1 a ij Yi’ < 
0 , where A = [a t j\ satisfies the comparison equation (13a). 

Proof We note that whenever V/(x,-( t)) = yf, for some i, and 
Vk(x k (r)) < Y k ^kYi, for some t > 0, we have 

Vi (x ,-)\ t=T < auY + {xk{Y)) < 0 • 

i.e. the (piecewise continuous) trajectories can never cross the 
boundaries defined as {x C S>\ V) (x,-) = yf Vf} . □ 

Propositions 1 and 2 can be used to replace the network-level 
conditions (11c) and (lid), respectively, by their equivalent 
decentralized, albeit more conservative, conditions to facilitate 
design of distributed control policies u, j Vi that satisfy 

Vi: W?(fi(xi)+u t j(xi)+gi{x)) < £ UjjVjiXj) Mxc2>, (14a) 

r'eAf 

f aij > 0 VjG,Af\{i}, 

subject to: < Y,jeJ< a ij an d (14b) 

HjeJi a ij 7 /* — 0 • 

4 In other words, a strictly diagonally-dominant matrix with negative diagonal 
entries is Hurwitz. 



Note that, a t j = 0 Vj (^Mf. Using the Positivstellensatz theorem 
(Theorem 2), with ki = (yf— V<(x,j) Vi, and VW=Qi, we can cast 
(14) into a set of SOS feasibility problems, for each i, 

-W?{fi+Ut,i+gi) + Yj ( a ij V j- a ij(y^- V j)) ( 15a ) 

je-A 

a5b) 

and e£ [°]- ( 15c ) 

where u t j G R[xj]”', VjGMf, (15d) 

a ,7 G R[0], and cijj G E[0] Vj G Mf\{i} . (15e) 

Here R[0] denotes scalar variables, £[0] denotes non-negative 
scalar variables and x, were defined in (6). 


The set of SOS conditions (15) defines the control u t l G R[x/]"‘ 
as an n,-vector of polynomials in x,, of a chosen degree. But 
further restrictions can be imposed on the control design. In 
this work, we consider bounded control signals of the form 


Vi: 

where, 


u t ,i,k(xi) | <U uk Vt > 0, Vk G {1,2,.. 

f Wt,i ~ Mt,i,2 ••• i 

[Uit> 0 V*g{1,2,. 


«/} (16a) 
(16b) 


For the uncontrolled states, we set the corresponding control 
bounds to zero. Further, by declaring these bounds as design 
variables the control problem can be formulated as a minimiza¬ 
tion of the maximal control efforts as, 

Vi: minimize Y}"' t / ( -a (17a) 

u t,i,ki a ij lk 


s.t., conditions (15), (17b) 

Uuk-u^k-G^itf-Vi) G £[*,■], V*G{1,...,/»,•}, (17c) 

U uk + GE[x,-], Vk G{1,...,*,■}, (17d) 

, - J = 0 , for the uncontrolled states, . 

w ere ’ i,k j G E[0 ], for the controlled states. e 

and o' j J /low G £[x;] V* G {1,2,. (17f) 

Given a choice of the degree of the control polynomials and an 
initial condition, (17) can be solved to find optimal, distributed, 
and exponentially stabilizing control policies. Algorithm 1 out¬ 
lines the major steps in the proposed control design procedure. 


It should be noted that for the subsystems that do not need to 
apply control the solution of the optimization (17) would result 
in Uj k = OVk G {1,.. 


Algorithm 1 Distributed Stabilizing Control Design 

procedure One-time Computation 
for each subsystem i G {1,2,..., m} do 
Compute the LF V)(x,j based on (7) 
Communicate Vj to neighbors M; Vj G ./tf\{;} 
Receive and store the LFs Vj(xj) Vj£jVi 
end for 

end procedure 

procedure Real-time Computation 
for each subsystem f G {1,2,..., m} do 
Compute initial level-set yf = V)(x,(0)) 
Communicate yf to neighbors 5?j Vj£^ij\{i} 
Receive y? Vj G {*} from neighbors 
Solve (17) for the optimal control input Utj(xj) 
end for 

end procedure 


Remark Often in practical scenarios, the control bounds need 
to be strictly imposed due to physical considerations, in which 
case the degree of the control polynomials can be varied to find 
feasible control policies. 

5. EXAMPLE 

We consider a network of nine Van der Pol ‘oscillators’ (see 
Van der Pol (1926)), with parameters of each oscillator chosen 
to make them individually (exponentially) stable (without the 
control). Each Van der Pol oscillator is treated as an individual 
subsystem, with the interconnections as shown below, 

: {1,2,5,9} : {2,1,3} M 3 : {3,2,8} 

M}: {4,6,7} MS: {5,1,6} {6,4,5} (18) 

M 7 : {7,4, 8 ,9} M6:{8,3,7} M$ : {9,1,7}. 

Each subsystem Vi G {1,2,... ,9} has two state variables, 
Xi = [x,m x,‘ 2 ] • The subsystem dynamics, under the presence 
of the neighbor interactions and control input, is given by 

: X[ 1 =X ()2 , 

Xi,2 = «;Xi,2 (l -Xi, 1 + Uf,i,2 + XjJ Y, PikX k ,2 ■ 

keJ'hli} 

where the subsystem parameters a, G [—2, —1] Vi and the inter¬ 
action parameters j3 j k G [—0.8,0.8] Vi , Vk G Mf\ {/}, are chosen 
randomly. Note that, we have considered u t .i.\ = 0 Vt Vi, i.e. the 
state variables x ,-1 Vi are not (directly) controlled. 

The goal is to apply the Algorithm 1 to compute distributed 
optimal controllers u t ,- 2 (jc,- i,x,- 2 ) Vi that guarantee exponential 
stabilization of the network of Van der Pol systems. 

5.1 Pre-Computation of Lyapunov Functions 

At first, we compute polynomial Lyapunov functions for the 
isolated (interaction free) and control-free subsystems 

Vi: x/,i =x, o, (20a) 

Xi ,2 = cCi Xj, 2 (1 -X 1 0 -X/, 1 , (20b) 

using the expanding interior algorithm (Section 3.1). As an 
example, we show a quadratic Lyapunov function and the 
associated estimate of the ROA of the interaction-free and 
control-free subsystem M 9 , 

= {(.X9,i,x 9 , 2 )\ V 9 < 1}, (21a) 

where, V9 = 0.595x9 .1 + 0 . 227 x 9 , 1 x 9,2 + 0.520x9 2 ■ (21b) 

Fig. 1 shows a comparison of the estimated ROA using the 
quadratic LF in (21), another estimate using a quartic LF 
and the ‘true’ ROA computed numerically by simulating the 
isolated and free dynamics. Clearly, the estimate improves with 
higher order LFs. However, for computational ease, the rest of 
the analysis will be based on quadratic LFs. 

Note that these LFs are computed only once for the network, 
and stored to be used for real-time control design. 

5.2 Controller Design: Test Case 

Figure 2 shows the evolution of the system state variables 
(belonging to subsystems Mj, M 3 , M 4 , M 5 , M 7 and Mg) and 
the subsystem LFs, starting from an unstable initial condition. 
In particular, the state variables belonging to the subsystems 
M 3 , M 7 and Ms ‘escape’ to infinity while other subsystems 
remain reasonably bounded, over the shown time window. 






(b) 


Fig. 2. System states (selected) and Lyapunov functions starting from an unstable initial condition, without any control. 



Fig. 1. Estimated ROA for isolated and free subsystem ,Zg. 

Algorithm 1 is used to compute distributed stabilizing linear 
controllers (with U t \ = 0 Vi), satisfying (15). Table 1 lists the 
results, while the trajectories after applying control are shown 
in Fig. 3. Interestingly, even though J ?3 was unbounded without 
control (Fig. 2), the algorithm finds that there is actually no need 
for control in .Zj provided its neighbors .Zj and ,9% remain 
bounded by their initial level-sets (Fig. 3). On the other hand, 
,9\ and ,9\ apply control, although they were bounded for over 
t € [0,2) without control (Fig. 2). 

The distributed control design is, however, conservative. For 
example, the maximum row-sum of the resulting comparison 
matrix (with control) is only marginally negative (Table 1), 
while its maximum eigenvalue actually turns out to be —0.06. 


Table 1. Distributed Control Results 


i 


E,=i a ‘j 


Ui, 2 

u t,i,2 i x i,l , x i,l) 

i 

0.33 

-0.000 

-0.029 

0.17 

-0.22.ri,! — 0.10x1,2 

2 

0.08 

-0.140 

-0.000 

0.00 

— 

3 

0.46 

-0.016 

-0.005 

0.00 

— 

4 

0.25 

-0.000 

-0.011 

0.12 

—0.18x4,1 —0.08x4,2 

5 

0.18 

-0.048 

-0.007 

0.00 

— 

6 

0.08 

-0.101 

-0.000 

0.00 

— 

7 

0.26 

-0.094 

-0.000 

0.65 

—0.47x7,1 — 0.89x7,2 

8 

0.94 

-0.000 

-0.109 

3.88 

—0.56x8,1 — 2.89x8,2 

9 

0.29 

-0.020 

-0.005 

0.00 

— 


6 . CONCLUSION 

The paper presents a distributed control strategy in which 
agents (subsystems) coordinate with their immediate neighbors 
to compute optimal local control strategies that exponentially 
stabilize the full nonlinear network. The proposed algorithm 
can be easily scalable to very large-scale, sparse, interconnected 
systems. Future work will explore ways to make the algorithm 
less conservative. One such way is to use a hierarchical two- 
level multi-agent control scheme, where the agents exchange 
some minimal information with a higher-level central agent. 
The central agent can perform minimal computations such as 
checking if the comparison matrix is Hurwitz (instead of the 
diagonally-dominant condition). Higher order polynomials for 
the subsystem Lyapunov functions could be used for potentially 
improved control design. It would be interesting to apply the 
proposed algorithm on some real-world system models, such as 
a network preserving power system network. 
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